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Abstract 

We establish a long time soliton asymptotics for a nonlinear system of wave equation coupled to a charged 
particle. The coupled system has a six dimensional manifold of soliton solutions. We show that in the large 
time approximation, any solution, with an initial state close to the solitary manifold, is a sum of a soliton 
and a dispersive wave which is a solution to the free wave equation. It is assumed that the charge density 
satisfies Wiener condition which is a version of Fermi Golden Rule, and that the momenta of the charge 
distribution vanish up to the fourth order. The proof is based on a development of the general strategy 
introduced by Buslaev and Perelman: symplectic projection in Hilbert space onto the solitary manifold, 
modulation equations for the parameters of the projection, and decay of the transversal component. 
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1 Introduction 



Our paper concerns an old mathematical problem of nonlinear field-particle interaction. A charged particle 
radiates a field which acts on the particle etc. This self-action is probably responsible for some crucial features 
of the process: asymptotically uniform motion of the particle, increment of the particle's mass etc. (see [38], 
Part I). The problem has many different appearances: for a classical particle coupled to a scalar or Maxwell 
field, for coupled Maxwell-Schrodinger or Maxwell-Dirac equations, for the corresponding second-quantized 
equations etc. 

One of the main goals of mathematical investigation of this problem is studying soliton-type long time 
asymptotics and asymptotic stability of soliton solutions. The first results in this direction have been discovered 
for KdV equation and other complete integrable equations. For KdV equation, any solution with sufficiently 
smooth and rapidly decaying initial data converges to a finite sum of solitons moving to the right, and a 
dispersive wave moving to the left. A complete survey and proofs can be found in [12^ 114]. 

For non-integrable equations, the long time convergence of the solution to a soliton part and dispersive 
wave was obtained first by Soffer and Weinstein in the context of ?7(l)-invariant Schrodinger equation with 
potential, |33 1 1341 [35] . The extension to translation invariant equations was obtained by Buslaev and Perelman 
[H m El E] for ID Schrodinger equation, and by Miller, Pego and Weinstein for ID modified KdV and RLW 
equations, [281 ISHl |30]. Later the results were developed by Bambusi and Cuccagna [H [HI [9l [10], and Martel 
and Merle [27]. 

In [H [5l [6] [7] the long time convergence is obtained for translation invariant and ?7(l)-invariant nonlinear 
Schrodinger equation: for any finite-energy solution ip{x,t) with initial data close to a soliton iPvq{x — v^t — 
oo)e*^°*, the following asymptotics hold: 

^p{x,t) =^py^{x-v±t-a±)e'^^^ + Wo{t)ip±+r±{x,t), t ^ ±oo. (1.1) 

Here the first term of the right hand side is a soliton with parameters v±, a±, u}± close to vq, oq, wq, the 
function W(){t)'ip± is a dispersive wave which is a solution to the free Schrodinger equation, and the remainder 
r±(x,t) converges to zero in the global L^-norm. 

In the present paper we consider a scalar real-valued wave field ip{x) in IR^, coupled to a relativistic particle 
with position q and momentum p, governed by 

■ip{x,t) = 7r{x,t), 7r{x,t) = Ail){x,t) — p{x — q{t)), x E IR^, 

. (1-2) 

q{t) = p{t)/^/l+W, p{t) = / V(x, t) Vp{x - q{t))dx. 

This is a Hamilton system with the Hamilton functional 

'H{'ip,iT,q,p) = ^ j (|7r(x)|^ + |V^(x)|2)(ix + j i){x)p{x - q)dx + ^/l+p^ . (1.3) 

The first two equations for the fields are equivalent to the wave equation with the source p{x — q). We have 
set the mechanical mass of the particle and the speed of wave propagation equal to one. The case of the point 
particle corresponds to p{x) = 6{x) and then the interaction term in the Hamiltonian is simply 'ip{q). However, 
in this case the Hamiltonian is unbounded from below which leads to the ill-posedness of the problem, that 
is also known as ultraviolet divergence. Therefore we smooth the coupling by the function p{x) following the 
"extended electron" strategy proposed by M. Abraham for the Maxwell field. In analogy to the Maxwell- 
Lorentz equations we call p the "charge distribution" . Finally, the form of the last two equations in ()1.2p is 
determined by the choice of the relativistic kinetic energy 
Let us write the system (jl.2p as 

Yit) = FiY{t)), t G IR, (1.4) 
where Y{t) := {'ip{x,t),TT(x,t),q{t),p{t)). The system (11. 2p is translation-invariant and admits soliton solutions 

Ya^v{t) = {tl^vix - vt - a),TTy{x - vt - a),vt + a,p^), p^ = v/y^l - v'^ (1.5) 
for all a,v ^ IR^ with \v\ < 1 (see ()2.8p to ()2.1ip below). The states Sa,v ■= ^a,t)(0) form the solitary manifold 



S := {Sa,v : a,veM^,\v\ < 1}. 



(1.6) 



Our main result (announced in [20]) is the soliton asymptotics of type (jl.ip . 



(V'(x,t),7r(x,t)) ~ {tpv^ix - v±t - a±),TTv^{x - v±t - a±)) + VFo(t)*±, t ±00, (1.7) 

for solutions to ()1.2p with initial data close to the solitary manifold S. Here Wo{t) is the dynamical group of 
the free wave equation, ^'j- are the corresponding asymptotic scattering states, and the remainder converges to 
zero in the global energy norm, i.e. in the norm of the Sobolev space see Section 2. For the 

particle's trajectory we prove that 

q{t) ^ v±, q{t) '-^ v±t + a±, t — )• ±00. (1.8) 

The results are established under the following conditions on the charge distribution: p is a real-valued 
function of the Sobolev class H^(]R^), compactly supported, and spherically symmetric, i.e. 

p,Vp,VVp£ L'^{M^), p{x) = for \x\ > Rp, p{x) = pi{\x\). (1.9) 

We require that all "nonzero modes" of the wave field are coupled to the particle. This is formalized by the 
Wiener condition 

p{k) = {2T^)-^'^ j e'^^ p{x)dx ^ {) for all A: G IR^ \ |q| _ (I JO) 

It is the nonlinear Fermi Golden Rule for our model: the coupling term p{x — q) is not orthogonal to the 
eigenfunctions e*'^^ of the continuous spectrum of the linear part of the equation (cf. [21 [6l [lOl [HI [32l [36]). 
Note that the Wiener condition is close to the linear version of the FGR | 3T1 pp. 67-68]. As we will see, the 
Wiener condition (ll.lOp is very essential for our asymptotic analysis. Generic examples of the coupling function 
p satisfying ()1.9p and (jl.lOp are given in |22j . In particular, the Wiener condition allows us to identify the 
discrete spectral subspace corresponding to the spectral point A = of the linearized system. Thus, we do not 
impose any implicit spectral conditions for the linearized system. Further, we will assume that p{k) has a fifth 
order zero at the point = 0, i.e. 

p '■"•'(O) = for all multiindeces a with \a\ < 4. (l-H) 

Equivalently, the following momenta vanish: 

j x"p{x)dx = 0, \a\ < 4, (1.12) 

in particular, the total charge of the particle is zero (neutrality condition). It is easy to obtain an example of 
p satisfying both (jl.lOp and (|l.lip . Indeed, take a p2 satisfying the Wiener condition (|1.10p . Then p = A^p2 
satisfies both (fLTnl) and (fLTTD . 

The system ()1.2p describes the charged particle interacting with its "own" scalar field. The asymptotics 
(jl.7p - (jl.8p mean asymtotic stability of uniform motion, i.e. "the law of inertia". The stability is caused by "ra- 
diative damping" , i.e. radiation of energy to infinity appearing analytically as a local energy decay for solutions 
to the linearized equation provided by the Wiener condition (II. 9p . The radiative damping was suggested first 
by M.Abraham in 1905 in the context of Classical Electrodynamics, [Ij. However, the asymptotics p.7p - p.8p 
are not proved yet for the Maxwell-Lorentz equations though close results are establihsed in |24j and |25j . 

One could also expect asymptotics ()1.7p for small perturbations of the solitons for the relativistic nonlinear 
wave equations and for the coupled nonlinear Maxwell-Dirac equations whose solitons were constructed in [3] 
and |13j respectively. Our result is a model of this situation though the relativistic case is still open problem. 

Let us briefly comment on earlier results. In the case of weak coupling, i.e. ||/o||l2 <C 1, scattering behavior 
of type (I1.7P for the system (II. 2p is established in \l7] for all finite energy solutions. In [15l [Ml [IS] the result 
was extended to the cases of Klein-Gordon field. Maxwell field, and spinning charge subject to Maxwell field 
respectively. The results under the Wiener condition are not established yet. 

The system ()1.2p under the Wiener condition was considered in |22t [23| 121] . In |22] the convergence to 
stationary states is proved for all finite energy solutions: in particular, the relaxation of accelleration g — )• as 
t — >• holds. In |23] the soliton-type asymptotics for the fields is established for all finite energy solutions. In 
|21] the effective dynamics under the slowly varying potential is constructed for solutions sufficiently close to 
the solitary manifold. However, the asymptotics p.7p and (II. 8p are not established in [22 | [23 | [21]. 



A long time asymptotics of type (jl.7p appears also in nonlinear wave equations, like the KDV [12\ I14j and 
the U (l)-invariant nonlinear Schrodinger equation [H O El [71 [29l [301 [Ml [Ml IM]- In these equations there are 
no particle degrees of freedom and the solitons (jl.5p correspond to the solitary wave solutions travelling at 
constant velocity. 

Our approach relies on and further develops the general strategy introduced in the cited papers in the context 
of the C/(l)-invariant Schrodinger equation. The approach uses i) symplectic projection of the dynamics in the 
Hilbert phase space onto the symplectic orthogonal directions to the solitary manifold to kill the runaway 
secular solutions, ii) the modulation equations for the motion along the solitary manifold, and iii) freezing 
of the dynamics in the nonauthonomous linearized equation. See more details in Introduction |19j where the 
general strategy has been developed for the case of the Klein-Gordon equation. The case of wave equation 
(jl.2p differs significally from the Klein-Gordon case because of i) slow Coulombic decay of the solitons, and ii) 
presence of the embedded eigenvalue in the continuous spectrum of the linearized equation (see the comments 
below) . 

Developing the general strategy for equations (II. 2p . we obtain our main result in Sections 3-6 and 11-15 
of the paper. The main novelty in our case is thorough establishing the appropriate decay of the linearized 
dynamics in Sections 7-10 and Appendices A, B: 

I. We do not postulate any spectral properties of the linearized equation, calculating all the properties from 
the Wiener condition (ll.lOp . Namely, we show that i) the full zero spectral space of the linearized equation 
is spanned by the tangent vectors, and moreover, ii) there are no others (nonzero) discrete eigenvalues (see 
Lemmas 19.51 19.61 and Proposition 18. ip . 

II. Using these spectral properties, we prove that the linearized equation is stable in the symplectic orthogonal 
complement to the tangent space Ts spanned by the tangent vectors da^Sa^v and dy.Sa,v, j = 1, 2, 3. We exactly 
calculate in Lemma 19.61 the corresponding symplectic orthogonality conditions for initial data of the linearized 
dynamics. 

III. One of main peculiarities of the wave equations (II. 2p is the presence of embedded eigenvalue A = in 
the continuous spectrum cic = IR of the linearized equation. This situation never happens in all previous works 
on the asymptotic stability of the solitary waves for the Schrodinger and Klein-Gordon equations. Thus, the 
symplectic orthogonality condition is imposed now at the interior point of the continuous spectrum in contrast 
to all previous works in the field. Respectively, the integrand at this point in the spectral representation of the 
solution is not smooth even if the symplectic orthogonality condition holds. Hence, the integration by parts in 
this spectral representation as in the case of the Schrodinger and Klein-Gordon equation, is impossible. For the 
proof of the decay in this new situation, we transform the spectral representation in the proofs of Propositions 
19.21 and [9.41 and develop new more subtle technique of convolutions. 

Our paper is organized as follows. In Section 2, we formulate the main result. In Section 3, we introduce 
the symplectic projection onto the solitary manifold. The linearized equation is defined and studied in Sections 
4-5. In Section 6, we split the dynamics in two components: along the solitary manifold, and in transversal 
directions, and we justify the slow motion of the longitudinal component. Sections 7-10 concern the time decay 
of the transversal component in the linearized dynamics. The time decay of the transversal component in the 
nonlinear dynamics is established in Sections 11-14. In Section 15 we prove the main result. In Appendices A, 
B we collect routine calculations. 

Acknowledgements The authors thank V. Buslaev for detailed and numerous lectures on his results, 
H. Spohn and E. Kopylova for fruitful discussions. 

2 Main Results 

2.1 Existence of Dynamics 

To formulate our results precisely, we need some definitions. We introduce a suitable phase space for the Cauchy 
problem corresponding to ()1.2p . Let be the real Hilbert space L^(IR^) with the scalar product (•,•) and 
the corresponding norm [ • [, and let be the completion of the real space C^(IR^) with the norm [V'(/'(2;)|. 
Equivalently, using Sobolev embedding theorem, = {^(x) G L^(IR^) : |V^(x)| G L^}- Let us introduce 
the weighted Sobolev spaces and with the norms \ip\a = 1(1 + |3;|)"^[ and = + |VV'[a 

respectively. 



Definition 2.1 i) £ is the Hilbert space © © IR'^ © IR^ with the finite norm 

II = IVV'I + |vr| + |g| + IpI for Y = {'il:,TT,q,p) . 

a) £a is the space © © IR^ © IR'^ with the norm 

\\Y\\a = II Y\\e^ = ||^/;||i,„ + I^U+i + \q\ + \p\. (2.1) 

Hi) T is the space © of fields F = {tp,TT) with the finite norm 

ll^lb = |V^| + |7r|. (2.2) 

iv) Toe is the space © L'^_^_i with the norm 

II^IU = ll^b. = llV'lli.a + lvrU+i (2.3) 

Note that is not contained in and for instance = oo if the neutrahty condition is not imposed, see 
(j2.1ip below. However, £ is the space of finite energy states (i.e. %{Y) < oo for Y G £) due to the fohowing 
estimates which are vahd for an arbitrary smooth i/>(x) vanishing at infinity 

"i// dxdy^^^^ = ^{p,A-'p)<^\V^P\' + {^{x),pix-q))< | V^l^ - A" V). (2.4) 

The Hamilton functional H is continuous in the space £, and the lower bound in (|2.4p implies that the energy 
(II. 3p is bounded from below. Note that the latter is not true if p is delta-function. 
We consider the Cauchy problem for the Hamilton system ()1.2p which we write as 

Y{t) = F{Y{t)), tGlR; Y{0)=Yo. (2.5) 

Here Y{t) = {ip{t),TT{t),q{t),p(t)), Yq = {iPot'^o,Qo,Po), and all derivatives are understood in the sense of 
distributions. 

Proposition 2.2 [22] Let i flTgl) hold. Then 

(i) For every Yq G £, the Cauchy problem \2. 5]) has a unique solution Y{t) G C(IR, f). 

(ii) For every t G IR, the map U{t) : Yq i— )• Y{t) is continuous in £. 
(Hi) The energy is conserved, i.e. 

n{Y{t))=n{YQ), tGlR, (2.6) 

and the velocity is hounded: 

\q{t)\<v<l, tGlR, (2.7) 

with some v which depends on Yq. 

2.2 Solitary Manifold and Main Result 

Let us compute the solitons (jl.Sp . The substitution to (jl.2p gives the following stationary equations, 

V 

Then the first two equations imply 

For \v\ < 1 the equation (j2.9p defines a unique function G i/^(IR^). If v is given and \v\ < 1, then p^ can be 
found from the third equation of (12. Sp . Further, functions p and ipy are even due to (|1.9p . Thus, V^v is odd 



Pv 



-v-\7TTy{y) = Aijy{y) - p{y) 

= -J ^My)p{y)dy 



(2.8) 



-A + {vVf)ij,{y) = -p{y), yGlR^ 



(2.9) 



and the last equation of ()2.8p holds. Hence, the soliton solution (jl.Sp exists and is defined uniquely for any 
couple {a,v) with \v\ < 1. 

The soliton can be computed by the Fourier transform tjj{k) := (27r)~'^/^ J e^^^^{x)dx: 



In the coordinate space 



V'^(x) = --3- [ -r— P^y^'^ y TT^{x) = -V -Vllj^ix), Pv=JV. (2.11) 

4vr7 h{y - x)ii + {y - x)±\ 

Here we set 7 = — v"^ and x = xy + x±, where x\\\\v and xx-Lu for x £ IR^. From the condition (jl.l2p it 

follows that 

ipv{y) ~ \y\~^, T^viy) ~ \y\~'^ as |y| 00 

and thus, 

V'^ e F^, a < 9/2; 7r„ G ^ a < 11/2. (2.12) 

Definition 2.3 A soliton state is S{a) := {ipv{x — b),T:^{x — b),b,pv), where a := {b,v) with b G IR^ and 
\v\ < 1. 

By (I2.12P for the soliton states we have 

S{a)G£a, a<| (2.13) 
Obviously, the soliton solution (II. 5p admits the representation Ya^vit) = S{a{t)), where 

a{t) = {b{t),v{t)) = {vt + a, v). (2.14) 



Definition 2.4 The solitary manifold is the set S := {S{a) : b G IR^, \v\ < 1}. 

The main result of our paper is the following theorem. 

Theorem 2.5 Let /il.9\) . the Wiener condition /il.lO\) . and the condition hold. Let < 6 < 1/2, set 

j3 := 4 + (5. Consider the solution Y{t) to the Cauchy problem h2. 5|) with the initial state Yq which is sufficiently 
close to the solitary manifold: 

YQ = S{ao) + Zo, d/3 := llZoll/3 < 1. (2.15) 
Then the asymptotics hold for t — t- ±00, 

m = v± + 0{\tr'''), qit)=v±t + a± + Oi\t\~^'), (2.16) 

(V'(x, t), 7r(x, t)) = (V'^i (x - v±t - a±), 7r^± (x - v±t - a±)) + Wo{t)^± + r±(x, t) (2.17) 

with 

\\r±{t)\\^ = 0{\t\-'). (2.18) 
It suffices to prove the asymptotics (I2.16p . (12.17P for t — )• +00 since the system (jl.2p is time reversible. 



3 Symplectic Projection 

3.1 Symplectic Structure and Hamilton Form 
The system (II. 2p reads as the Hamilton system 



Y = JVniY), J :-- 



( 





1 









-1 
























V 








-/3 


/ 



(3.1) 



where DTi is the Frechet derivative of the Hamilton functional ()1.3p . I3 is the 3x3 identity matrix. Let 
us identify the tangent space to £, at every point, with £. Consider the bilinear form defined on £ by 



where 



dipi^x) A dn^x) dx + dq A dp, i.e. 



n{Yi,Y2) = {Y^,JY2), Yi,Y2e£, 



(3.2) 



{Yi,Y2) := (-01,^2) + (7ri,7r2) + qiq2 + P1P2 
and stands for the scalar product / ■ipi{x)il^2{x)dx or its different extensions. 



Definition 3.1 i) Yi \ Y2 means that Yi is symplectic orthogonal to Y2, i.e. Q(Yi,Y2) = 0. 

a) A projection operator P : £ ^ £ is called symplectic orthogonal i/Yi f Y2 for Yi G KerP and I2 £ Im-P- 

3.2 Symplectic Projection onto Solitary Manifold 

From now on we suppose that the condition (jl.lip is satisfied. Let us consider the tangent space Ts(^a)S to the 
manifold 5 at a point S{a), where a = {b,v). The vectors tj := dajS{a), where d^^ ■= dh^ and 8^^^^ := d^^ 
with j = 1, 2, 3, form a basis in 7^5. In detail, 



Tj = Tj{v) := db^S{a) = ( -djtp^iy) , -djTry{y) , ej , ) 
T3+j = T3+j{v) := dy^S{a) = ( d^^^^{y) , dy^TTy(y) , , d^^p^ ) 



j = 1,2,3, 



(3.3) 



where y := x — b is the "coordinate in the moving frame", ei = (1,0,0) etc. 
By (j2.12p for the tangent vectors we have 



9 

Tj{v) e£a, a< -; j = 1, 



,6. 



(3.4) 



Lemma 3.2 The matrix with the elements n{Ti{v),Tj{v)) is non- degenerate for \v\ < 1. 



The proof is made by a straightforward computation, see [19], Lemma 3.2 for the case m = (note that the 
left hand side of the identity [191 (A. 8)] is well defined by the condition (jl.lip ). 

Now we show that in a small neighborhood of the soliton manifold S a "symplectic orthogonal projection" 
onto S is well defined. Let us introduce the translations Ta : {ip{x),Tr{x),q,p) 1— )• {■ip{x — a),7r(x — a),q + a,p), 
a G IR^. Note that the manifold S is invariant with respect to the translations. 

Definition 3.3 Let us denote by v{Y) := p/\/l where p G IR'^ is the last component of the vector Y. 
Lemma 3.4 Let lil.9\) hold, a > —9/2 and v < 1. Then 

i) there exists a neighborhood Oa{S) of S in £a and a map H : 00,(5) — )• S such that H is uniformly continuous 
on Oa{S) n {Y G £a ■ v{Y) < v} in the metric of £a, 



ny = y for YeS, and Y-S\TsS, where S = IiY. 
ii) Oa{S) is invariant with respect to the translations Ta, and 

IlToY = TJTY, for Y G Oa{S) and a G IR^. 
Hi) For any v < \ there exists a v <1 such that \v{YlY)\ < v when \v{Y)\ < v. 



(3.5) 
(3.6) 



iv) For any ij < 1 there exists an ra{v) > s.t. S{a) + Z G Oa{S) if \v{S{a))\ < v and \\Z\\a < ra{v). 



The proof is similar to that of [19', Lemma 3.4]. 

We will call 11 the symplectic orthogonal projection onto S. 
Corollary 3.5 [191 Corollary 3.5]. The condition \2.15\) implies that Yq = S + Zq where S = TIYq, and 

\\Zoh < 1. (3.7) 

4 Linearization on the Solitary Manifold 

Let us consider a solution to the system ()1.2p . and split it as the sum 

Y{t) = S{a{t)) + Z{t), (4.1) 

where a{t) = (6(t),z;(t)) G IR^ x {| v\ < 1} is an arbitrary smooth function of t G IR. In detail, denote 
Y = {il;,7T,q,p) and Z = {'^,U,Q,P). Then (jUD means that 

i^{x,t) = V,(t)(x-6(t)) + 1'(x-6(t),t), q{t) = b{t)+Q{t) 
7T{x,t) = Tr,^t)ix-bit))+U{x-bit),t), Pit) = p,(t)+P(t) 

Let us substitute (|4.2p to (jl.2p . and linearize the equations in Z. Later we will choose 5((7(t)) = IlY{t), 
i.e. Z(t) is symplectic orthogonal to Ts(a{t))'S- However, this orthogonality condition is not needed for the 
formal process of linearization. The orthogonality condition will be important in Section 6, where we derive 
"modulation equations" for the parameters ait). 

Let us proceed to linearization. Setting y = x — bit) which is the coordinate in the moving frame, we obtain 
from (j42]) and ([O]) that 

V ■ V^^^(t) iy)-b- VV'^(t) iy) + ^iy,t)-b- V^iy, t) = 7r„(j) (y) + U{y, t) 

V ■ "^vT^vit) iy)-b- V7r^(t) iy) +Uiy,t)-b- VU{y, t) 
AiP^^t){y) + A^{y,t) - p{y - Q) 

Pv + P 



TT 



b + Q 



Vl + (p, + P)2 



(4.3) 



p = v-VvPv{t)+P = -{^i'<Pv{t)iy) + ^iy.t)),piy-Q)) 

The equations are linear in ^ and IT, hence it remains to extract linear terms in Q and P. Within this section, 
we estimate the remainders in the norms of the space £a with an arbitrary a > 0. 

First note that piy - Q) = piy) - Q ■ Vp(y) - iV2(Q), where -iV2(Q) = piv - Q) - piv) + Q ■ ^piv); for 
N2iQ) the bound holds, 

\N2iQ)\a < CiQ)Q^ (4.4) 
uniformly in \Q\ < Q for any fixed Q. Second, the Taylor expansion gives 

^"^^ =v + uiP- viv ■ P)) + Nsiv, P), 

Vl + ip^ + P)2 

where u = := (1 +pl)^^^'^ = ^/\ — v^, and 

\N-iiv,P)\<Civ)P'' (4.5) 

uniformly in v with \v\ < v < 1. Using the equations (|2.8p . we obtain from (j4.3p the following equations for 
the components of the vector Z{t): 

^(y, t) = n(y, t)+b- V^iy, t) + {b - v) ■ VMv) " v ' "^vMv) 



Uiy,t) = A^iy,t)+b-VU{y,t)+Q-Vp{y) + ib-v)-VTr^iy)-v ■V^TTyiy) + N2 
Qit) =Uy{l3-v(g)v)P + iv-b) + N3 



(4.6) 



Pit) 



= (^iy,t),Vpiy)) + {VMy), Q- Vp(y)) - v ■ VyPy + NJv, Z) 



where iV4(^^, 2') = {Vipy, N2{Q)) + {V'^ ,Q-V p) + {V'^ , N2{Q)) . Clearly, Ni{v,Z) satisfies the following estimate 



\N^{v, Z)\ < Cp{p, V, Q) + ll^-H.^IQl] , 

uniformly in v^Q with \v\ < v and \Q\ < Q. We can write the equations (I4.6P as 

Z{t) = A{t)Z{t)+T{t) + N{t), tGlR. 



(4.7) 



(4. 



Here the operator A{t) depends on a{t) = {b{t),v{t)). We will use the parameters v = v{t) and w := b{t). 
Then A{t) can be written in the form 



A{t) 



n 

Q 
\P J 



An' 



/^\ 


/ w-V 


1 


\ 


/^\ 


n 


A 


W-V V/9- 




n 


Q 


■~ 


OB,, 




Q 


\P J 




(vv^,-Vp) 


/ 


\P J 



(4.9) 



where B^ = v^uil^ — v ® v) . Furthermore, T{t) and N{t) in (j4.8p stand for 



T{t) = 



I {W - V) ■ Vlpy - V ■ VvlPv \ 
{W — V) ■ VVTt, — V ■ VyTTy 

V — w 

-V ■ VvPv 



, N{t) = N{a,Z) 



( \ 

N2{Z) 

N^{v,Z) 
\ N^{v, Z) ) 



(4.10) 



where v = v{t), w = w{t), a = a{t) = {b{t),v{t)), and Z = Z{t). Since \Q\ < for any q, the estimates 

(fO]) with Q = r_„({;), (j43]) and (|47l) imply the following 



(4.11) 



Lemma 4.1 For any a > 

\\N{a,Z)\U<C{v)\\Z\\l 
uniformly in a,Z with \\Z\\^a ^ f-a{v) and \v\ < v. 



Remarks 4.2 i) The term A{t)Z{t) in the right hand side of the equation (j4.8p is linear in Z{t), and N{t) is 
a high order term in Z{t). On the other hand, T(t) is a zero order term which does not vanish at Z(t) = 
since S{a{t)) generally is not a soliton solution if (|2.14p does not hold (though S{a{t)) belongs to the solitary 
manifold). 

ii) Formulas ()3.3p and ()4.10p imply: 



T{t) = - ^[(w - v)lTi + VlTi+s] 



(4.12) 



1=1 



and hence T(t) £ 75(o-{t))'5, t G IR. This fact suggests an unstable character of the nonlinear dynamics along 
the solitary manifold. 



5 The Linearized Equation 

Here we collect some Hamiltonian and spectral properties of the operator (14. 9p . The statements of the section 
are particular cases of those in [19], Section 5 for m = 0. First, we consider the linear equation 



X{t) = Ay^y,X{t), t G IR 
with an arbitrary fixed v such that \v\ < 1, and w G IR^. Let us define the space 

= H'^iM^) e H^OR^) e ir2 e irI 



(5.1) 



Lemma 5.1 i) For any v, \v\ < 1 and w G IR^ the equation /i5. 1]) formally can be written as the Hamilton 
system (cf. 

X{t) = JDn,o,u,{Xit)), t £ M, (5.2) 
where DTiy^^ is the Frechet derivative of the Hamilton functional 



n,,UX) = \j [|n|2 + |V^|2]dy + j TlwV^dy + j p{y)Q- 



V^dy 

+ lp.B,P-l(g.VV.(y),g-Vp(y)), X = {^,Xl,Q,P)e£. (5.3) 

a) Energy conservation law holds for the solutions X{t) G C(Si,£), 

nv,w{X{t)) = const, t £ M. (5.4) 

Hi) The skew-symmetry relation holds, 

n{A,,^Xi,X2) = -n{Xi,A,^^X2), x^gs, X2££+. (5.5) 

Lemma 5.2 The operator A^^^ acts on the tangent vectors Tj{v) to the solitary manifold as follows, 

Ay^^[Tj{v)] = {w -V)- VTj{v), A„,^[Tj+3(u)] = {w -v) ■ VTj+3{v) + Tj{v), j = 1, 2, 3. (5.6) 

We will apply Lemmas 15.11 and 15.21 mainly to the operator A^^^ corresponding to w = v. In that case (|5.5p 
reads 

AvA-^ji^)] = 0, A^A'^j+siv)] = Tj{v), j = 1,2,3. (5.7) 
Moreover, the linearized equation acquires the additional essential feature. 

Lemma 5.3 Let us assume that w = v and \v\ < 1. Then the Hamilton functional \5. 3|) reads, see [19' (5.14)] 

n^A^) = ^ / (|n + ^ • v^-p + |Ay2^ - A-^/2g . vp|2^ dx + ]^p- b^p > o. (5.8) 



Here A^, is the operator defined by h2.9\) . 

Remark 5.4 Lemma [5.31 together with energy conservation (j5.4p imply the analyticity of the resolvent (.A-^ — 
A)^^ for ReA > 0, see below. 

Remark 5.5 For a soliton solution of the system (jl.2p we have 6 = u, u = 0, and hence T{t) = 0. Thus, the 
equation (j5.ip is the linearization of the system (jl.2p on a soliton solution. In fact, we do not linearize (jl.2p on 
a soliton solution, but on a trajectory S{a{t)) with a{t) being nonlinear in t. We will show later that T(t) is 
quadratic in Z{t) if we choose S{a{t)) to be the symplectic orthogonal projection of Y{t). Then (|5.ip is again 
the linearization of (11.21). 



6 Symplectic Decomposition of the Dynamics 

Here we decompose the dynamics in two components: along the manifold S and in transversal directions. The 
equation (|4.8p is obtained without any assumption on a{t) in ()4.ip . We are going to choose S{a{t)) := ny(t), 
but then we need to know that 

Y{t) G Oa{S), t G IR, (6.1) 

with some (Da{S) defined in Lemma |3.4[ It is true for t = and a < (3 := 4 + 5 hy (j3.7p . if > in (|2.15p is 
sufficiently small. Then ^(^(O)) = ny(0) and Z(0) = Y{0) - S{a{0)) are well defined. 

We set a = —f3 and will prove below that (|6.ip holds if djs is sufficiently small. First, the a priori estimate 
(12:711 together with LemmaE31iii) imply that nY{t) = S{a{t)) with a{t) = {b{t),v{t)), and 

\v{t)\<d<l, t£M (6.2) 

if Y{t) G 0^(5). Denote by ra{v) the positive number from Lemma 13.41 iv) which corresponds to the chosen 
a = —f3. Then S[a) + Z G 00,(5) if o" = {h,v) with < v and \\Z\\a < ra{v). Note that (j3.7p implies 
||Z(0)||a < ra{v) if d/3 is sufficiently smah. Therefore, S{a{t)) = nY{t) and Z{t) = Y{t) - S{a{t)) are well 
defined for t > so small that ||Z(i)||Q, < ra{v). This is formalized by the following standard definition. 



Definition 6.1 is the "exit time", 

= sup{t > : < r_^(i)), < s < t}, Z{s) = Y{s) - S{a{s)). (6.3) 

One of our main goals is to prove that = oo if is sufficiently small. This would follow if we show that 

\\Z{t)\\_p <r_(^{v)/2, 0<t<t,. (6.4) 
Now N{t) in (j4.8p satisfies, by (|4.1ip with a = /3, the following estimate: 

\\Nit)\\0<C{v)\\Z{t)\\lf,, 0<t<t,. (6.5) 

6.1 Longitudinal Dynamics: Modulation Equations 

From now on we fix the decomposition Y{t) = S{a{t)) + Z{t) for < t < t^, by setting S{a{t)) = IlY{t) which 
is equivalent to the symplectic orthogonality condition of type (|3.5p . 

Zit)\Tsi^it))S, 0<t<t,. (6.6) 



This allows us to simplify drastically the asymptotic analysis of the dynamical equations (j4.8p for the transversal 
component Z{t). As the first step, we derive the longitudinal dynamics, i.e. the modulation equations for the 
parameters o'{t). Let us derive a system of ordinary differential equations for the vector cr(t). For this purpose, 
let us write (16.61) in the form 



Q(Z(t),Tj(t)) = 0, i = 1,...,6, 0<t<t*, (6.7) 
where the vectors Tj(t) = Tj{a{t)) span the tangent space Ts(a{t))'S- Note that a{t) = {b{t),v{t)), where 

\v{t)\<v<l, 0<t<t^, (6.8) 

by Lemma 13.41 iii). It would be convenient for us to use some other parameters {c,v) instead of o" = {b,v), 

where c{t) = h{t) — I v{T)dT and 
Jo 

c{t) = b{t) -v{t) =w{t) -v{t), 0<t<t^. (6.9) 
We do not need an explicit form of the equations for (c, v) but the following statement. 

Lemma 6.2 LetY{t) be a solution to the Cauchy problem 112. 5\) . and lll.ll]) . ((^. ( [g. 7| ) hold. Then {c{t),v{t)) 
satisfies the equation 

m 



y -AA(a(t),Z(t)), 0<t<t„ (6.10) 
where 

M{a,Z) = 0{\\Z\t^) (6.11) 

uniformly in a ^ {{b^v) : \v\ < v}. 

Proof We differentiate (j6.7p in t and take the equation (j4.8p into account. Then (see details of computation 
in [19], Lemma 6.2) we obtain, in the vector form [191 (6.18)]: 

= n{v) (^".^ +Moia,Z) J ^ +Moia,Z), Mojia, Z) = niN,T,). (6.12) 

Here the matrix il.{v) has the matrix elements il.{Ti,Tj) and hence is invertible by Lemma [3.21 The 6x6 matrix 
Mo{a,Z) has the matrix elements ~ ||Z||_^ and hence we can resolve the equation ()6.12p with respect to 
{c,v). Then (j6.1ip follows from Lemma |4 . 1 1 with a = /3, since A/q = OdlZH^o). □ 



6.2 Decay for the Transversal Dynamics 



In Section 15 we will show that our main Theorem 12.51 can be derived from the following time decay of the 
transversal component Z(t): 

Proposition 6.3 Let all conditions of Theorem \2. 51 hold. Then = oo, and 



i^(t)ii-^< :;^':'7; , t>o. (6.13) 



We will derive (|6.13p in Sections 8-14 from our equation (j4.8p for the transversal component Z{t). This 
equation can be specified using Lemma W% Namely, by (j4.12p and ()6.9p 

3 

T{t) = - ^[cm + U/Tz+s]. 

(=1 

Note that the norm ||T(t)||^ is well-defined by p.4p . Then Lemma 16.21 implies 

\\T{t)\\p<C{v)\\Z{t)\tp, 0<t<t,. (6.14) 

Thus, in (|4.8p we should combine the terms T(t) and N{t) and obtain 

Z{t) = A{t)Z{t) + N{t), 0<t<t^, (6.15) 

where A{t) = ^^(t),^(t), and N{t) := T{t) + N{t). From (i6?T4]l and (|63]l we obtain that 

\\N{t)y<C{v)\\Z{t)f_p, 0<t<t,. (6.16) 

In the remaining part of our paper we will analyze mainly the basic equation (I6.15P to establish the decay 
()6.13p . We are going to derive the decay using the bound ()6.16p and the orthogonality condition ()6.6p . 

Let us comment on two main difficulties in proving (I6.13p . The difficulties are common for the problems 
studied in O |8]. First, the linear part of the equation is non- autonomous, hence we cannot apply directly 
known methods of scattering theory. Similarly to the approach of [H |8] , we reduce the problem to the analysis 
of the frozen linear equation, 

X{t) = AiX{t), t e IR, (6.17) 

where Ai is the operator Ay^^y^ defined in ()4.9p with vi = v{ti) and a fixed ti G [0,t*). Then we estimate 
the error by the method of majorants. Let us note that recently some methods of freezing were developed by 
Cuccagna and Mizumachi, [10, Jlj. 

Second, even for the frozen equation (j6.17p . the decay of type (|6.13p for all solutions does not hold without 
the orthogonality condition of type (16. 6p . Namely, by ()5.7p the equation (I6.17P admits the secular solutions 

3 3 

X{t) = ^CjT,{vi) + Y,Djh{vi)t + Tj+s(.vi)] (6.18) 
1 1 

which arise also by differentiation of the soliton (jl.Sp in the parameters a and vi in the moving coordinate 
y = X — vit. Hence, we have to take into account the orthogonality condition (j6.6p in order to avoid the secular 
solutions. For this purpose we will apply the corresponding symplectic orthogonal projection which kills the 
"runaway solutions" (I6.18p . 

Remark 6.4 The solution (I6.18P lies in the tangent space Ts(ai)<S with ai = {bi,vi) (for an arbitrary bi G IR) 
that suggests an unstable character of the nonlinear dynamics along the solitary manifold (cf. Remark 14. 21 ii)). 

Definition 6.5 i) Denote by Uy, \v\ < 1, the symplectic orthogonal projection of £ onto the tangent space 
Ts{a)S, and Py = 1- Ily. 

ii) Denote by Zy = Py£ the space symplectic orthogonal to Ts(cj)S with a = {b,v) (for an arbitrary b £ M,). 



Note that by the Hnearity, 

U^Z = Y,^dy)rj{vMTiiv),Z), Ze£, (6.19) 

with some smooth coefficients Yljii^v^ . Hence, the projector n^, in the variable y = x — b, does not depend 
on 6, and this explains the choice of the subindex in 11^, and P^. Now we have the symplectic orthogonal 
decomposition 

£ = Tsia)S + Zv, <y = {h,v), (6.20) 
and the symplectic orthogonality (16. 6p can be written in the following equivalent forms, 

n,(t)Z(t) = 0, P,(t)Z{t) = Z{t), 0<t<U. (6.21) 

Remark 6.6 The tangent space Ts(a)<S is invariant under the operator A^^^ by Lemma 15.21 hence the space 
is also invariant by (jS.Sp : A^^^Z £ for sufficiently smooth Z £ Z^,. 

The following proposition is one of main ingredients for proving (I6.13p . Let us consider the Cauchy problem 
for the equation (|6.17p with A = A^^^ for a fixed v, \v\ < 1. Recall that the /3 = 4 + (5, 0<5< 1/2. 

Proposition 6.7 Let /il.9\) . the Wiener condition /il.lO\) . and the condition hold, \vi\ < v < 1, and 

Xq £ £. Then 

i) The equation 1^6.11 ), with Ai = A^-^^y^, admits the unique solution c^'^^Xq := X{t) G Cb(Si,£) with the initial 



condition X{0) = Xq. 

a) If Xq £ Zy^ n the following decay holds, 

||e^^%||-2-5 < JYJ^^Uoh^ t £ M. (6.22) 

Part i) follows by standard arguments using the positivity (j5.8p of the Hamilton functional. Part ii) will be 
proved in Sections iTlflOl developing general strategy [I9l [20] . Namely, the equation (j6.17p is a system of four 
equations involving field components, ^ and H as well as vector components, Q and P. We apply Fourier- 
Laplace transform, and express the field components in terms of the vector components. Then we obtain a 
closed system for the vector components alone and prove their decay. Finally, for the field components we come 
to a wave equation with a right hand side which has the established decay. This implies the corresponding 
decay for the field components. 

7 Solving the Linearized Equation 

Here we start the proof of Proposition 16.71 solving the linearized equations ()6.17p . Let us apply the Laplace 
transform 

poo 

AX = X{X) = / e-^^X{t)dt, ReA > (7.1) 

to (|6.17p . The integral converges in since ||X(t)||£: is bounded by Proposition 16.71 i). The analytic- 
ity of X{\) and Paley-Wiener arguments should provide the existence of a i5-valued distribution X{t) = 
{^{t),U{t),Q{t),P{t)), t £ IR, with a support in [0,oo). Formally, 

A^^X = X{t) = — [ e''^^X{iu; + 0)duj, t £ M. (7.2) 
In 



2tt 



To prove the decay (j6.22p . we have to study the smoothness of X(ia; + 0) at G IR. After the Laplace transform 
the equation (j6.17p becomes 

\X{\) = AiX{X) + Xq, ReA>0. 
In detail (for simplicity we write v instead of vi in this section), 

n(y) + V ■ V§(y) - A§(y) = -§o(y) 



A^r(y) + V ■ Vn(y) + Q ■ Vp{y) - Xfi{y) = -Ilo{y) 

ByP -\Q = -Qo 
-(V*(y), p(y)) + (VMy), Q ■ Vp(2/)) - AP = -Po 



y £ IR3. (7.3) 



Let us consider the first two equations. In Fourier space they become 

U{k) -ivk^ik) - Xif{k) = -^o{k) 



-k^^{k) - {ivk + X)tl{k) = -Uo{k) + iQkp{k) 



/c G (7.4) 



This imphes 



^ = i{{ikv + X)4>o + Uo-ikQp), (7.5) 
n = i(-A;2^o + [ikv + A)no - i{ikv + \)kQp), (7.6) 



where 



D = D{X) = k'^ + {ikv + Xf. (7.7) 

From now on we use the system of coordinates in x-space in which v = (|f|,0,0), hence vk = \v\ki. Substitute 
(j7.5p to the 4-th equation of ()7.3p and obtain 

j ^-^{{ikv + A)^o + no - ikQp)f)dk + j kip^kQjdk - XP = -Pq. 

Since ipy = —p/{k'^ — (kv)'^), we come to 

{K - H{X))Q + XP = Po + $(A). 

Here 

^>(A) = $(^&o,no)(A):=i I i{{ikv + X)4>o + Uo)Jdk = i{ ^'^"^ + + "° , kp), (7.8) 
and H(X) are 3 x 3-matrices with the matrix elements 

fhkjimidk f hk,\p{krdk 

The matrix K is diagonal and positive definite since p{k) is spherically symmetric and not identically zero 
by (jl.lOp . The matrix H is well defined for Re A > since the denominator does not vanish. The matrix H 
is diagonal similarly to K. Indeed, if i 7^ j, then at least one of these indexes is not equal to one, and the 
integrand in (|7.9p is odd with respect to the corresponding variable by (jl.Op . Finally the 3-rd and the 4-th 
equations of (|7.3p become 



M(^) ( ? ) = ( P, (A) ) ■ *^(^> = ( -nx) 'xil ) ■ "^'^ = - 

Remark 7.1 Note that 

$(A) = ^^0, no)(A) = A{W\t){^o, Ho), Vp), (7.11) 

where W^{t) is the first component of the dynamical group W{t) defined below by (jlO.Sp and A is the Laplace 
transform (17.11). 



Let us proceed to ^-representation. We invert the matrix of the system (j7.4p and obtain 



-k^ -{ivk + X) J i\ ' J ' i ]^ f.2 _(^iyk + X) 



Taking the inverse Fourier transform, we obtain the corresponding fundamental solution 

Gx{y)=('"^~^ ^, , )9x{y), (7.12) 



-A vV - X ^ 

where g\{y) is the unique tempered fundamental solution of the determinant 

D = D{X) = -A + {-vV + Xf. (7.13) 

Thus, 

9xiy) = F^Xy -^—^-^—^ = F-^y -^—^-^-^——^ , y e IRl (7.14) 
Note that the denominator does not vanish for Re A > 0. This implies 



Lemma 7.2 The operator G\ with the integral kernel G\{y — y') is continuous operator H^{M?) © L'^{R^) 
F2(IR3) h^(r3^ for ReA > 0. 



From (|7.5p and (|7.6p we obtain the convolution representation 

^ = _(^.V-A)c/A*^o + 5A*no + (5A*Vp)-Q 

n = Agx*^o-iv'^-^)9x*^o-iv-V-X)igx*Vp) 



(7.15) 



Let us compute gx{y) explicitly. First consider the case v = 0. The fundamental solution of the operator 
-A + A2 is 

gx{y) = ^Yj. (7.16) 

Thus, in the case v = we have 

-A -1 \ e-^\y'y'\ 



-A J 47r|y — y'\ 

For general v = (|f|,0,0) with \v\ < 1 the denominator in (j7.14p . which is the Fourier symbol of D, reads 

D{k) = k^ + {i\v\ki + Xf = {l-v^)kl + kj + kl + 2i\v\kiX + X^ 

i\v\X 

where 



(1 - v'){k^ + + ki + kl + x^ (7.17) 



9,2X2 X2 

x2 = f^ + A2 = -^. (7.18) 



Therefore, setting 7 := l/\/l — t)2, we have 

x = 7A. (7.19) 



Return to x-space: 



1 



D = -(yi+jxiY + ^i:=-f\v\X (7.20) 

T 



Define yi := 72/1 and Vi := d/dyi. Then 

D = -(Vi + j<i)2- V|- V^ + x^ (7.21) 

Thus, its fundamental solution is 



e 



-K\y\-3<iyi 



9\{y) = — r~Fi — ' y := (72/1,2/2,2/3)- (7.22) 

By (I7T9I1 . (I7:20]l we obtain 

< Re xi < Re X, Re A > 0. (7.23) 
Let us state the result which we have got above. 

Lemma 7.3 i) The operator D = D{X) is invertible in L^(]R^) for ReA > and its fundamental solution 



7.2S^ decays exponentially in y. 



a) The formulas {7.2!^ and ( 7.19), ( 7. 2(J^ imply that the distribution g\{-) admits an analytic continuation 
in the papameter X from the domain Re A > to the entire complex plane (D. 

Lemma 7.4 The matrix function M (X) (respectively, M~^{X)) admits an analytic (respectively meromorphic) 
continuation from the domain Re A > to the entire complex plane (D. 



Proof From the first equation of (j7.15p and the last equation of ()7.3p it fohows that 

H,,{X) = {g^*d,p,d,p) (7.24) 

and thus, by ([722]), 

/(J 
_ djp{x)djp{y)dxdy < oo. 
l-*^ y\ 

By dZSI), (fTlOD this implies 

sup |F(A)| < oo. (7.25) 

RcA>0 

The analytic continuation of M(A) exists by the expressions (|7.24|) and Lemma 17.31 ii). since the function p{x) 
is compactly supported by ()1.9p . The inverse matrix is then meromorphic since it exists for large Re A. The 
latter follows from (fTlOD since H{X) ^ 0, ReA oo, by ([73]). □ 

8 Regularity in Continuous Spectrum 

Proposition 8.1 The matrix M^^{iuj) is analytic in a; € IR\ {0}. 

Proof It suffices to prove that the limit matrix M{iuj) := M{iuj + 0) is invertible for a; 7^ 0, a; G IR if /) satisfies 
the Wiener condition (jl.lOp . and |?;| < 1. Let > 0. One has 

det M{iuj) = det ( '"^i]. , 7^"" ^ = -(a;2 + u^fi){uj^ + uff, w e IR, (8.1) 
y —t'izuj) tujl-i J 

where F[ioj) := F{iuj + 0), /i := Fii{iijj), and / := F22{iuj) = i<33(ia;). The invertibility of M{iuj) follows from 
(jS.ip by the following lemma, whose proof is based on the Sokhotsky-Plemelj formula, see [i39i. Chapter VII, 
formula (58)]. 

Lemma 8.2 // il.lO\) holds, then for a; G IR the imaginary part of the matrix - — -F{iuj) is negative definite, 

\uj\ 

i.e. - — 7lmFjj{iu}) <0, j = 1,2,3. 

\uj\ 

Proof Since F{ioj) = H{iuj + Q) — K, where the matrix K is real, we will consider only the matrix H(iuj + 0). 
For e > we have 



f k'Ap{k)\^dk 

+ ^) = y J- = 1. 2. 3. (8.2) 



Consider the denominator 

b{iuj + e,k) = k"^ - {\v\ki + uj - ief . 
D{iuj, A;) = on the ellipsoid T^^ if |a;| > 0, where 

T^ = {k:{uk,-^-^r + ki + kl = ^}, 



here u = \J\ — . From the Sokhotsky-Plemelj formula for C^-functions it follows that 



Imi7,-,(ia; + 0) = -— TT / — dS, (8.3) 

1^1 JT^\VD{iuj,k)\ 

where dS is the element of the surface area. This immediately implies the statement of the Lemma since the 
integrand in (j8.3p is positive by the Wiener condition (jl.lOp . This completes the proofs of the lemma and the 
Proposition 18.11 □ 

Remark 8.3 The proof of Lemma 18.21 is the unique point in the paper where the Wiener condition is indis- 
pensable. 



9 Time Decay of the Vector Components 

Here we prove the decay (|6.22p for the vector components Q{t) and P{t) of the solution c'^^^Xq. Formula (|7.10p 
expresses the Laplace transforms Q{X), P{X). Hence, the components are given by the integral 

V P{t) J 2ttJ '\ Po + Htuj) J ^ ' 

Let us recall that in Proposition 16.71 ii) we assume that 

XQeZ^r\£p, /3 = 4 + (5, 0<5<l/2. (9.2) 



Theorem 9.1 The functions Q{t), P{t) are continuous for t > 0, and 

+ (TTSyLll^oll/3, t>o. (9.3) 

Proof Note that the Proposition 18.11 alone is not sufficient for the proof of the convergence and decay of the 
integral ()9.ip . We need an additional information about the regularity of the matrix M~^{iuj) at its singular 
point oj = 0, and some bounds at |a;| — >• oo. 

Let us split the integral ()9.ip in two terms using the partition of unity Cii^) + C2('^) = 1, G H: 

Pit) )-2^J' (^^^"^ + I P{^u) J - I Pi(t) J + I P2{t) ) ' ^'-'^ 

where the functions Ck{^) G C°°(]R) are supported by 

suppCi C {w G IR : < r + 1}, suppC2 C {w S IR : > r}, (9.5) 

where r is introduced below in Lemma [931 We prove the decay (19. 3p for (Qi, -Pi) and {Q2, P2) in Propositions 
19.41 and 19.21 respectively. 



Proposition 9.2 The functions Q2(t), P2{t) are continuous fort > 0, and 

\Q2{t)\ + \P2it)\ < (Y^^^]Lll^oll/5- (9-6) 

Proof First we study the asymptotic behavior of Af~-^(A) at infinity. Let us recall that M^^(X) was originally 
defined for Re A > 0, but it admits a meromorphic continuation to (D (see Lemma [73]). 

Lemma 9.3 There exist a matrix Rq and a matrix-function Ri(uj), such that 

M-\iu;) = — + Riiu), |a;| > r > 0, a; G IR, (9.7) 

CO 

where, for every k = 0,1, 2, 

|a^i?i(a;)| < j%, |u;|>r>0, ujGM, (9.8) 

r is sufficiently large. 

Proof The statement follows from the expilicit formulas ()A.4p to (jA.SP for the inverse matrix M~^(iuj) and 

from the bound (|7.25p . □ 



Q2(.t) ^ _ ^ / e-*C2(^)M-i(i^) 



Further, (j9.ip imphes that 

Pq \ ^(lUj) 



P2{t) ; 27r 



m[ h (9.9) 



Po J V / 

where fsee (|7.2p ) 

s{t) := A"^ [C2{oo)M-\iu)] 

and 

/(t) := A-ici>(ic^) = {w\t){^oM,Vp), (9.10) 
since ^> is given by (j7.8p and (|7.1ip . Note that s{t) is continuous for t > 0, and 

\s{t)\ = 0{t-^), t ^ oo, y N > 0. (9.11) 



by (|9.7p - (|9.8p . On the other hand, recah that, under the conditions of Proposition 16.71 and Theorem I2.5| 
{^0, Hq) G J-jS with /3 = 4 + 6 where 5 > 0. Then we obtain that 

m\ < ^,%'L \\^oh (9.12) 

by Lemma [122] below with a = /3. Hence, dMI) fohows from ([Ml) by (IOT]l and IMTTIh . □ 

Now let us prove the decay for Qi{t) and Pi{t). In this case the proof will rely substantially on the symplectic 
orthogonality conditions. Namely, (19. 2p implies that 

0(Xo,r,) = 0, j = 1...6. (9.13) 



Proposition 9.4 The functions Qi{t) and Pi{t) are continuous for t > 0, and 



cip,v)\\X oy 



\Qiit)\ + iPi(t)i < :r, .1+7 > * > 0- (9-14) 



Proof First, we obtain the formulas for the Fourier transforms Q{iuj) and P{iuj) 
Lemma 9.5 The matrix M^^{iu) can be represented as follows, 



Liu) := M-\u) = ^ , (9.15) 



where Cij{uj), i,j = 1,2 are smooth diagonal 3 x 3-matrices, Cij{uj) G C°°(— r — l;r + 1), and 

Cii = iCi2B-\ (9.16) 
For proof see Appendix A. Now (j7.10p implies that the vector components are given by 

Q{iuj) = -Cii{u)Qo + \ci2{^){Po + ^{iuo)), (9.17) 

P{iu) = C2i{uj)Qo + -£22(a;)(Po + <^iito)), (9.18) 

CO 

Next we calculate the symplectic orthogonality conditions (j9.13p . 
Lemma 9.6 The symplectic orthogonality conditions jy.lS]) read 

Po + $(0) = and B-iQo + ^'(0) = 0. (9.19) 



For proof see Appendix B. 

Now we can prove Proposition 19.41 
Step i) Let us prove (j9.14p for Pi{t) relying on the representation ()9.18p . Namely, (j9.4p and ()9.18p imply 

Pi{t) = A-iCi(^)^2i(a;)Qo + A-\iico)C22{co)^^^^^ = PHt) + (t). 

CO 

The first term -P{(t) decays like Ct~°°||Xo||/3 by Lemma 19.51 The second term admits the convolution repre- 
sentation Pi{t) = A~^Ci£22 * g(t), where 

git) .= A-i^I±£M_ 
Now we use the symplectic orthogonality conditions (j9.19p and obtain 

= A-' = (.,.20) 

OO 

Therefore, P['{t) decays like Ct^^^+S) ||Xo||^ for t > 0, since by (I9J2D 

\g{t)\<Cip,v)il + t)-^^+'^\\Xoh, t>0. (9.21) 

Step ii) Now let us prove (|9.14p for Qi{t). By (j9.17p . (|9.16p . and the symplectic orthogonality conditions (|9.19p . 

Qiiuj) = zP 'Qq h = \iB^ ^Qo H = 

CO \ 00 J LO \ UJ J 

[iPv Qo+9{cO))=Cl2 = ^^12 , 

CO CO CO 

since iB^^Qq + ^(0) = by the symplectic orthogonality conditions (|9.19p . because ^(0) = i<I>'(0). Thus, 
Qiit) = A-iCi(w)£i2 * h{t) by ([921), where 

t 

n{t) := A = I / g[s)ds, 

00 

similarly to i^lH^i . This integral decays hke Ct-^-^||Xo||,3 for t > by dOT]) . hence (l9lH) for Qi{t) is proved. 
The proof of Proposition 19.41 and Theorem 19.11 is complete. □ 

10 Time Decay of Fields 

Here we construct the field components ^{x, t), n(x, t) of the solution X{t) and prove their decay corresponding 
to (j6.22p . Let us denote F{t) = {'^{■,t),I[{-,t)). We will construct the fields solving the first two equations of 
(j6.17p . where A is given by (j4.9p . These two equations have the form 

-(')Ka" v,.,,))- 

By Theorem 19.11 we know that Q{t) is continuous and 



cip,v)\\X o\y 



\Q{t)\< r/Vr ^ i>o. (10.2) 



Hence, the Proposition 16.71 is reduced now to the following 



Proposition 10.1 i) Let a function Q{t) £ C([0, oo); IR^), and Fq G F. Then the equation \10.1\) admits a 
unique solution F[t) G C[0,cxd; J^) with the initial condition -F(O) = Fq. 

a) If Xq = {Fq]Qq,Pq) G 8p and the decay il0.2\) holds, the corresponding fields also decay uniformly in v: 



for \v\ < V with any v G (0; 1). 

Proof The statement i) follows from the Duhamel representation 



(10.3) 



Fit) = Wit)Fo + 



[ W{t-s) 
Jo 





Vp-Q{s) 



ds 



t > 0, 



where W{t) is the dynamical group of the modified wave equation 



vV 1 



F{t). 



A vV 

The group W{t) can be expressed through the group Wo{t) of the wave equation 



1 
A 



Namely, the problem (|1U.6|) corresponds to pU.Sp . when v = 0, and it is easy to see that 

[W{t)F{0)]{x) = [Wo{t)F{0)]{x + vt), a; G ]R^ t G IR. 



(10.4) 



(10.5) 



(10.6) 



(10.7) 



The identity (jl0.7p implies the energy conservation law for the group W{t): for {^{■,t),Il{-,t)) = W{t)F{0) 
we have 

\U{x, t)-v V^'(x, t)\^ + I V^'(2;, t)p) dx = const, t G IR. 



(10. 



/ 



In particular, this gives by (j2.2|) 

\\W{t)F4:F<C{v)\\Fo\\T, iGH. 
This estimate and (jl0.4p imply the statement i). 
Let us proceed to the statement ii). 

Lemma 10.2 For v < \ and Fq G J-'a, a > I, the following decay holds, 



t > 0, 



(10.9) 



for the dynamical group W{t) corresponding to the modified wave equation 1^10. 5|] with \v\ < v. 



Proof For the case v = the proof is provided in [26]. For a nonzero v with |w| < {; the proof is similar, we 
provide it for convenience. 

We should estimate ||VF(t)Fo||_Q for large t > 0. Set e = (1 — v)/2. For an arbitrary sufficiently large t > 1 
let us split the initial function Fq in two terms, Fq = Fq^ + Fq ^ such that 



l|-^0,J|o + ll-^o'tlU ^ C'll-^olU) t 1, 
where C does not depend on t, and 

FQt{x) = 0, |x| > et, FQt{x) = 0, |x| <et-l. 



(10.10) 



(10.11) 



For an arbitrary / G H"^ and a > 1 one has < see [2^, formula (2.9). Now the estimate for 

W{t)F^^ follows by (fnTH]) . (fToTTl) and ffTOTO]) : 



\\W{t)F'\\_o. < C\\Wit)Fn\^ < lb < Ci(^)ll^aJa(l+N)"" < C^2(i))||Fo|U(l+N)-", t > 1. (10.12) 



It remains to estimate W{t)FQ^. First note that 

Wo{t)F^ t{x) = for l^l < (1 - e)t (10.13) 

by the strong Huygen's principle for the group Wo{t). The principle reads 

Woix-y,t)=0, \x-y\^t, (10.14) 

where Wo{z, t) is the integral (distribution) matrix kernel of the operator VFo(t). Further, from ()10.13p it follows 
that 

[Ty(t)Fo_J(x) =0 for |x| < et 

by (110. 7p and since \v\ < v = 1 — 2e. 

For an arbitrary / G such that f{x) = in the region {|x| < et}, one has < C'('u)t~""'"^||/||^i, see 

[26] . the proof of Proposition 2.1. Applying to / = W{t)FQ ^ we obtain by (|10.8p that, 

\\W{t)F^^^ < C{v)t-'^^'\\W{t)F^4^ < a-"+i||F^,,|b < Ct-'^^'\\F^,t\\T^ < Ct-"+i||Fo|b,. 

The proof is complete. □ 

Now the statement ii) of Proposition [Tim follows from the Lemma [10.2l and the Duhamel representation (|10.4p . 
Indeed, 

\\W{t)Fo\\-2-s < Ct-^-^WFoh+s < Ct-^^^\\Xo\\2+s < Ct-^-^\\Xo\\p 
by Lemma 110.21 with a = 2 + 6. Further, 

/ll(t),Vp-Q(s))||2+jds 







{l + it-s)) 



1+S 



t t 
, /■ \Q{s)\ds „ r ds 

- J {i + (t-s)y+^- " °"''y (i + (t_5))i+5(i+,)i+5 



by Lemma 110.21 with a = 2 + 6, regularity properties of p, and (jl0.2p . The last integral decays like (1 + t)~^~^ 
by a well known result on decay of a convolution. □ 
The proof of Proposition 16.71 is complete. □ 

11 Frozen Form of Transversal Dynamics 

In next four sections we will prove the transversal decay (j6.13p relying on the Proposition 16.71 First, let us fix 
an arbitrary ti E [0,t*), and rewrite the equation (j6.15p in a "frozen form" 

Z{t) = AiZ{t) + {A{t)-Ai)Z{t) + N{t), 0<t<t„ (11.1) 

where Ai = and 



A{t) - Ai 



/ [w{t)-v{ti)] • V \ 

[w{t)-v{ti)]-V 

B,^t)-B,^t,) 



\ (V(V'.(t)-^.(t0)rVp) / 

The next trick is important since it allows us to kill the "bad terms" [w{t)—v{ti)] ■ V in the operator A{t) — A\. 

Definition 11.1 Let us change the variables {y,t) i— )• {yi,t) = (y + di{t),t) where 

di{t) := [ {w{s) -v{ti))ds, 0<t<ti. (11.2) 
Jti 



Next define 

Zi{t) = {^i{yi,t),Ui{yi,t),Q{t),P{t)) := {^{y,t),U{y,t),Qit), Pit)) 

= {^{yi-di{t),t),U{yi-di{t),t),Q{t),P{t)). (11.3) 

Then we obtain the final form of the "frozen equation" for the transversal dynamics 

Zi{t) = AiZi{t) + Bi{t)Zi{t) + Ni{t), 0<t<ti, (11.4) 

where Ni{t) = N{t) from the basic equation ()6.15p expressed in terms of y = yi — di{t), and 



Bi{t) 



/ \ 



B^(t-^-B^(t^-) 

V (V(^,(t)-V',(t,)),-Vp) J 



At the end of this section, we will derive appropriate bounds for the "remainder terms" Bi{t)Zi(t) and Ni{t) 
in (jll.4p . First, note that we have by Lemma 16.21 



|5.(t) - < I / Hs) ■ V,i?,(,)ds| < C / ' \\Z{s)f_f,ds. (11.5) 

Jti Jt 

Similarly, 

KV(V„(i) - V.(to)'-Vp)| < ' \\Z{s)\\l^ds. (11.6) 
Let us recall the following well-known inequality: for any a G IR 

(l + |?/ + x|)" < (l + |y|)"(l + |x|)l"l, x,yeM^. (11.7) 



Lemma 11.2 [19' Lemma 7.2]. For {^,Il,Q,P) £ with any a G IR the following estimate holds: 

||(^'(yi-di),n(yi-di),Q,P)|U < ||(M/,n,Q,P)||„(l + |(ii|)H , 4 G Ml (11.8) 
Corollary 11.3 The following bound holds 

\\N,it)y<\\Ziit)\\l^il + \di{t)\)'^, 0<t<ti. (11.9) 
Indeed, applying the previous lemma twice, once for a = /3 and once for a = — /?, we obtain from (j6.16p that 
||A^i(t)||^ < (1 + \dm)^\\Nit,Zit))\\^ < (1 + \dmf\\Z\tfs < (1 + \d,{t)\?^\\Zi{t)\tfs ■ 

Corollary 11.4 The following bound holds 

\\Bi{t)Zi{t)\\^ < C||Zi(t)||_^ f\l + \di{T)\)'('\\Zi{T)\tf,dT , < t < tl. (11.10) 



For proof we apply Lemma 111.21 with a = —/3 to ()11.5p and ()11.6p and use the fact that Bi{t)Zi{t) depends 
only on the finite-dimensional components of Zi{t). 



12 Integral Inequality 

Recall that < S < 1/2. The equation ()11.4p can be written in the integral form: 

Zi{t) = e^i*Zi(0) + / e^i(*-')[5iZi(s) + Ni{s)]ds, 0<t<ti. (12.1) 
Jo 

We apply the symplectic orthogonal projection Pi := Pv{ti) to both sides, and get 

PiZi(t) = e^i*PiZi(0) + /*e^i(*-^)Pi[SiZi(s) + iVi(s)]ds. 

Jo 

We have used here that Pi commutes with the group e"^^* since the space Zi := Pi£ is invariant with respect 
to e"^^*, see Remark 16.61 Applying (I6.22p we obtain that 



|PlZi(t)||_2-5 < 



c 



Jo 



since the operator Pi is continuous in Ep. Hence, from (112. 2p and (|11.9p . (|11.10p we obtain that 



(12.2) 



\\PiZv 
+ C(di 



-2-S 



< 



c 



1 



\Zi{0)h 



||Zi(.)||„^ / ||Zi(r)rVr+||Zi(.)||l^ 



ds, 0<t< ti, 



(l + |t-s|)i+5 

where di := supo<t<tj |'^i(t)|- Since ||Zi(t)||-|-^ < C((ii)||Z(t)||-|-^ by Lemma [11.21 '^e can rewrite ()12.3p as 

C(di' 



(12.3) 



|PlZi(t)||_2-5 < 



+ C{di) 



(1 + tY+^ 
1 



11^(0)11/3 







(l + |t_s|)l+5 
Let us introduce the majorant 



\\Z{s)\\^^ I \\Z{r)\\l^dr+\\Z{s)\tf, 

3, te[o,Q. 



ds, 0<t< ti, 



m{t) := sup {l + sy+^\\Z{s 
se[o,t] 



To estimate di{t) by m{ti) we note that 



w{s) — v{ti) = w{s) — v{s) + v{s) — v{ti) = c(s) + I v{T)dT 
by (|6.9p . Hence, pi.2p . Lemma 16.21 and the definition (jl2.5p imply 



(12.4) 



(12.5) 



(12.6) 



\di{t)\ 



{w{s) - v{ti))ds\ < 
ti / 1 



c s 



\v{T)\dT I ds 



(1 + 5)2+2-5 



+ 



dT 



(1 + r) 



2+25 



ds < Cm2(ti), < t < ti. 



(12.7) 



We can replace in (jl2.4p the constants C{di) by C if m{ti) is bounded for ti > 0. In order to do this 
replacement, we reduce the exit time. Let us denote by e a fixed positive number which we will specify below. 



Definition 12.1 is the exit time 

< = sup{t G [0, t^) : m(s) < e, < s < t}. 
Now (fTTip and (TTT} imply that for ti < t'^ 

C 



(12.8) 



|PlZi(t)||_2-5 < 



-c 



{i + tY+^ 



11^(0)11/3 



(l + |i-5| 



\l+5 



\\Z{s)\\.p / \\Z{T)\\ipdT+\\Z{s)\\tp 



ds, 0<t< ti. 



(12.9) 



13 Symplectic Orthogonality 

Finally, we are going to change PiZi{t) by Z{t) in the left hand side of (jl2.9p . We will prove that it is possible 
using again that djj <^ 1 in ()2.15p and due to the following important bound: 

Lemma 13.1 For sufficiently small e > 0, we have for ti < t'^; 

\\Z{t)\U-5<C\\PiZi{t)\\_2-5, 0<t<tu (13.1) 

where C depends only on p and v. 

Proof The proof is based on the symplectic orthogonality (|6.2ip . i.e. 

n,(t)Z(t) = 0, tG[0,ti], (13.2) 

and on the fact that all the spaces Z{t) := P,u{£^£ are almost parallel for all t. 

Namely, we first note that \\Z{t)\\-2-5 < C(e)||Zi(t)||_2-5 by LemmaHOl since \di{t)\ < Ce^ for t < ti < < 
by p2.7p . Therefore, it suffices to prove that 

||^l(t)||-2-5 < 2||PiZi(i)||_2-5, 0<t<ti. (13.3) 

This estimate will follow from 

\\n,,Zi{t)\\_2-5<\\\Zi{t)\\.2-S, 0<t<ti, (13.4) 

since PiZi{t) = Zi{t) — n^^Zi(t), where vi = v{ti). To prove (|13.4p . we write ()13.2p as 

n,(j),iZi(t) = 0, tG[0,ti], (13.5) 

where II^^^) x^i(0 is n^(^)Z(t) expressed in terms of the variable yi = y + di{t). Hence, (jl3.4p follows from 
(jl3.5p if the difference II^j — n^(j) x is small uniformly in t, i.e. 

||(n,, -n,(i),i)^i (011-2-5 < ^11^1 (t)||-2-5, o<t<ti. (13.6) 

It remains to justify (jl3.6p for a sufficiently small e > 0. In order to prove the bound ()13.6p . we will need the 
formula (|6.19p and the following relation which follows from (j6.19p : 

n,(t),iZi(t) = ^n,,(t;(t))T,-i(?;(t))f](T,,i(t;(t)),Zi(0), (13.7) 

where Tj^i{v(t)) are the vectors Tj{v{t)) expressed in the variables yi. In detail (cf. (|3.3p ). 

TjAv) ■■= i-dj'ilj^{yi-di{t)),-djTTy{yi- di{t)),ej,0), = 123 (138) 
Tj+3,i(w) := {d^ipy{yi - di{t)),d-a^7r^{yi - di{t)),0,dy^p.a), ' ' ' 

where v = v{t). Thus, we have to estimate the difference of 

n„,Zi(t) = ^IVji{vi)Tj{vi,yi)9.{Ti{vi,yi),Zi{t,yi)) 

and 

n,(t),iZi(t) = ^%(t;(t))Tj(^;(t),2/i - di{t))^{n{v{t),yi - di{t)), Zi{t,yi)). 
The estimate is based on the following bounds. First, 

lUjiivit)) -Uji{v{ti))\ = \ J^\is) ■V^nji{v{s))ds\ <C J^' \v{s)\ds, 0<t<ti, (13.9) 

since \VyIlji{v{s))\ is uniformly bounded by (|6.8p . Second, 

\^{Ti{vi,yi) - Ti{v{t),yi - di{t)),Zi{t,yi))\ < \\Ti{vi,yi) - Ti{v{t),yi - di{t))\\2+s\\Zi{t,yi)\\_2-d- 



Further, since < Ce^ and Vr^ are smooth and sufficiently fast decaying at infinity functions, Lemma 

111.21 implies 

- Tj{vm\2+s < C\di{t)\ < Ce\ 0<t<ti (13.10) 
for all j = 1, 2, . . . , 6, where C depends only on 5 and v. Finally, 

Tj{v{t)) - Tj{v{ti)) = / V{s) ■ V^Tj{v{s))ds, 



and therefore 



\Tj{v{t)) - rj{viti))\\2+5 < C / \ijis)\ds, 0<t<ti. 



(13.11) 



At last, the bounds (flalHI) will follow from (f6l9]) . (fT37fll and (ll3.inp - (fT3:9]l if we establish that the integral in 
the right hand side of (|13.11|) and (|13.9p can be made as small as we please by choosing e > sufficiently small. 
Indeed, 

ds 



\v{s)\ds < Cm\ti] 



{l + s 



i2+25 



< Ce, 0<t<ti. 



The proof is complete. 



(13.12) 

□ 



14 Decay of Transversal Component 

Here we complete the proof of Proposition [6? 



Step i) We fix an e, < e < r_^('u) and t'^ = t'^,(e) for which Lemma 113.1 1 holds. Then the bound of type (jl2.9p 
holds with ||PiZi(t)||_2-(5 in the left hand side replaced by ||Z(t)||_/3 : 



C 



\\Z{t)\\.p < \\Z{t)\\_2-5 < C\\PiZ,{t)\U^s < 



\z 



+c 



(1 + 1* 



\l+<5 



\Z{s)\\_^ \\Z{r)f_^dT+\\Z{s)\\l^ 



ds, <t <ti 



(14.1) 



for ti < t'^. This implies an integral inequality for the majorant m{t) introduced by ()12.5p . Namely, multiplying 
both sides of ()14.ip by (1 + 1)^^^, and taking the supremum in t G [0, ti], we get 



m{h)<C\\Z{0)y + C sup / 



mis) 



+ 



1+5 



ti 



m?{T)dT 



2+2<5 



+ 



m?{s) 



(1 + ^) 



2+25 



ds 



for ti <t'^^. Taking into account that m{t) is a monotone increasing function, we get 

m(ti) < C\\Z{mp + C[m3(ti) + m\ti)\I{ti), h < C 

where 

1 dr 1 



(14.2) 



(l + t)^+^ 

I{tl) = sup / 

telctilJo (l + |i-s|)^+'^ 



+ 



ds < I < oo, ti>0. 



(l + s)l+5 (1+^)2+25 (l + 5)2+25_ 

Therefore, ()14.2p becomes 

m{h) < C||Z(0)||^ + CI[m\h) + m'^ih)], h < C (14.3) 
This inequality implies that m{ti) is bounded for ti < t'^, and moreover, 

m(ti) <Ci||Z(0)||;3, h<t:, (14.4) 
since m(0) = ||Z(0)||^ is sufficiently small by (|3.7p . 

Step ii) The constant Ci in the estimate (|14.4p does not depend on and t'^ by Lemma [13.11 We choose 

in ([235]) so small that ||^(0)||/3 < e/(2Ci). It is possible due to dST])- Then the estimate pTI]) implies that 



t'^ = t^ and therefore ()14.4p holds for all ti < t^. Then the bound ()12.7p holds for aU t < t^. Therefore, ([67 
also holds for all t <t^. Finally, this implies that = oo, hence also = oo and (|14.4p holds for all ti > if 
dp is small enough. 

The transversal decay (|6.13p is proved. □ 



15 Soliton Asymptotics 

Here we prove our main Theorem 12 . 51 reiving on the transversal decay (j6.13p . First we will prove the asymptotics 
(j2.16p for the vector components, and afterwards the asymptotics (j2.17p for the fields. 

15.1 Asymptotics for the vector components 

From (113]) we have q = b + Q, and from ([6J5|) . ([6J6|) with /3 = 4 + 5, and (gj]) it follows that Q = 
B,^,)P + 0{\\Z\\l^). Thus, 

q = b + Q = v{t) + c{t) + B,^t)P{t) + 0{\\Z\\lf,). (15.1) 
The equation (j6.10p and the estimates (j6.1ip . (|6.13p imply 

m\ + m\<f0i, t>o. (15.2) 

Therefore, c{t) = c+ + C'(t-(i+2'5)) and v{t) = v+ + 0(^^(1+2'^)), t oo. Since |P| < \\Z\\_fs, the estimate 
(|6T3]1 . and ([TO]) . (fTET]l imply that 

qit)=v+ + Oit-^'^). (15.3) 

Similarly, 

b{t) = c{t)+l v{s)ds = v+t + a+ + 0{t-'^^), (15.4) 
Jo 

hence the second part of (|2.16p follows: 

q{t) = b{t) + Q{t) = v+t + a+ + 0{t-^^), (15.5) 

since Q{t) = 0{t-^-^) by (f03D . 



15.2 Asymptotics for the fields 

We apply the approach developed in [T7], see also [151 llSl 113 EI]- For the field part of the solution, F{t) = 
{'ip{x,t),TT{x,t)) let us define the accompanying soliton field as 

-^v{t)(i) = (V'v(t)(3; - q{t)), Tr^(t){x - q{t))), 

where we define now v(t) = q{t), cf. (jl5.ip . Then for the difference Z(t) = F{t) — i^v(t)(0 '^^ obtain easily the 
equation [21], Eq. (2.5), 

Z(t) = ^Z(t)-v-VvF,(i)(t), A(V,vr) = (7r,A^). 
Then ^ 

Z{t) = Woit)ZiO) - [ Wo{t - s)[v{s) ■ VvF,(,)(s)]ds. (15.6) 
Jo 

Since ||(V'i>+,7r„+)(a^ - v+t - a+) - F^(f)(t)||j- = 0{t-^^) by (fTO]) and pX5|) . to obtain the asymptotics ([2T7P 
it suffices to prove that Z{t) = WQ{t)^+ + r+{t) with some ^+ € T and ||?'+(i)||j- = 0{t~^). This is equivalent 
to 

Wo{-t)Z{t) = ^+ + r'4t), (15.7) 

where ||r^!,_(t)||jr = 0{t~^) since Wo(*) is a unitary group in the Sobolev space F by the energy conservation 
for the free wave equation. Finally, (]15.7p holds since (115. 6p implies that 

Wo{-t)Z{t) = Z{0) + /* Wo{-s)R{s)ds, R{s) = v(s) • VvF,(,)(s), (15.8) 
Jo 

where the integral in the right hand side of (|15.8p converges in the Hilbert space J- with the rate 0{t^^). 
The latter holds since ||Wo(— s)-R('S)||.F = 0{s^^^^) by the unitarity of Wo{—s) and the decay rate ||i?(s)||jr = 



0{s ^ Let us prove this rate of decay. It suffices to prove that |v(s)| = 0{s ^ or equivalently \p{s) 
0{s~^~^). Substitute (j4.2p to the last equation of (|1.2p and obtain 



pit) = I [ij,^t){x-bit)) + ^{x-b{t),t)]Vpix-b{t)-Q{t))dx (15.9) 



The first integral in the right hand side is zero by the stationary equations (j2.8|) . The second integral is 
0{t~'^-^), since Q{t) = 0{t-^-^), and by the conditions ([L9]) on p. Finally, the third integral is C'(t"^"'^) by 
the estimate (|6.13p . The proof is complete. □ 

A Structure of the matrix M~^{iu) 

We prove Lemmas 19.31 and 19.51 Recall that for G IR 

M{iuj) = 

where 



F{iuj) iujis 





















\ 





^ ) 





Moj) 

B^= \ ly \ , F{iuj) = | /(w) 

f{u) 

Here = Fii{iu + 0), /(w) = ^22(^6^ + 0) = ^33(^0; + 0) with 



\pm f i/5|2fc2 



' fc2 + (A + ikivy J k'^ - (kiv) 
Fjj{X) are analytic functions in (D by Lemma EH Thus, 

F,,(A) = F,,(0) + Fj,.(0)A + ^^A2 + . . . 



Here Fjj{0) = by (lAT]) . Further, by fOl 

(fc2 + (A + iffci)2)2 



F;,(A) = -2 / ciA;fc^|p|2-^4±^^^^i-^ (A.2) 



and 



since the integrand function is odd in ki. Hence, we obtain Fjj{X) = X'^rj^X), where ^^(A) is analytic in C 
Note that rj{0) = F^'j{0)/2. By fO]) we have 



Fj^.(A) = -2jdkk]\p\ 



2| -12 - 3(A + it)A;i^2 



(fc2 + (A + zz;A:i)2)3 
and finally. 



Let us denote r{uj) = r2(w) = r3(a;). Then 



^^"^^='^ (^")=U2lM ^22H J' 



(A.4) 



where 



/ — idJ 



Lii{uj) 



V 










-iuj 











-iuj 



1 



—I 

1 - v^r\{uj) 











—i 





1 — vr(u) 








— i 
\ — vr 



by (|A.3P : we denote the last matrix LwiuS). Similarly, 



Li2(a;) 


















1 - v^riiuj) 










1 — vr(uj) 




(A.5) 



1 — vriLo) ) 

(A.6) 



we denote the last matrix £12 (w). Note that 



Further, 



21 



/ -/i(a^) 








\ 





so we put £21 = -^21- Finally, 






-/(^) 



( n(^) 

1 - i^^ri(a;) 










1 — i/r(ci;) 










(A.7) 



1 — vr(bj) j 



L22{U}) = Lii{uj) = -Cu{uj) 



(A.8) 



and thus, £22 1'^) = Note that the denominators of the matrix elements of each matrix Cu to £22 are 

nonzero at a; = 0, since ri(0) < and r(0) < by (IA.3p . For a; 7^ the denominators are nonzero by Lemma[821 
This completes the proof of Lemma l9.5i Finally, ()A.4p to (jA.SP imply Lemma 19.31 since rj{uj) = —Fjj{iuj)/uj'^, 
where Fjj{iu) are bounded functions by (I7.25p . 

B Symplectic orthogonality conditions 



Let us check that the symplectic orthogonality conditions (]9.13p with j = 1,2,3 read the first equation of 
dllS]). By dZS]), 

$(M/o,no)(0) = i{ikv^o + fto, -^), ^(0) = k^- {kvf. 
On the other hand, by (l9l^ with j = 1, 2, 3, and ([33]), (I2J0I1 . we have 



Z)(0) 



D(0)' 



= {kv^,,^^)-^{U,,^^) - (Po), = -<f,(v&o,no)(0) - (Po),. 
Now let us check that the conditions (j9.13p with j = 4, 5, 6 read the second equation of (j9.19p . By (|7.8 



$'(0) = ^{^,kp) - ^{2^kv'M^,kp) = + ^""^'^^^^ + '^^"^ 



'D(O) 



P'2(o) 



(A;2 - {kvff 



,kp), 



where the integral converges by the condition (jl.lip . On the other hand, by (|9.13p with j = 4,5,6, and (|3.3p . 
(12^01) . we have for j = 1, 2, 3 



since Qo ■ d^.p^ = Qo ■ ^ej = ^Qq ■ ej. 
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